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1, INTRODUCTION 

Rapid increasing and most requests have to tend for the applications of BLDC motor started from 
household applications tends to industrial and aerospace services. As all the BLDC motor applications are due 
to its advantages such as its high dynamic response, wide speed range, high efficiency of operation and 
noiseless [1, 2]. 

BLDC motor is similar to the synchronous motor with the permanent magnet. Electronic switches 
inverter 1s used for current commutation. The inverter is supplied from dc supply that the name came. The rotor 
position, detected either by hall sensor or any other method, is necessary for deriving the inverter to supply the 
three-phase trapezoidal stator windings. No need to model the damper winding and the rotor is lapsed because 
of high resistivity of the stainless-steel magnet [3]. To make synchronous between back emf's and stator 
currents, the last have to be quasi-square shape to produce high torque and steady speed. The six-step strategy 
is the known operation for stator that makes two phases on and the third is float [4-6]. 

Different speed control techniques are proposed and implemented to achieve an accurate torque rate 
and speed control behavior. The control used in most research depends on the PID controller in which its 
parameters are tuned either by trial and error, Ziegler-Nicholas, a fuzzy optimization technique, genetic or 
particle swarm optimization [5-7]. Advanced techniques to improve the operation of the PID controller are also 
interested such as employing fractional-order PID-Controller and using this controller via two order of freedom 
[8-10]. Some researchers have expanded the structure of Artificial intelligence to control the nonlinear system 
such as wavelet and wavelet neural network (Wave net) approaches [11-14]. 

Speed and current limiting strategy make a necessary role in modern control for BLDC motor drive. 
This paper proposes a pulse width modulation (PWM) system for BLDC motor to control the speed and limit 
the current during starting and operating conditions using two PID controllers. Each PID controller parameters 
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are tuned using the proposed harmony search optimization technique. The MATLAB software program 1s used 
to show the complete system processes. The BLDC motor performance in terms of speed control and 
modification of current are tested and the results are discussed. 


2. EQUIVALENT CIRCUIT AND MODELING OF BLDC MOTOR 

As the stator of synchronous motor, the stator of BLDC motor appears. It has a three-phase winding 
in the stator while the rotor consists of permanent magnets. A style of BLDC motor supplied from the three- 
phase inverter is shown in Figurel. The ultimate voltage frontier of the stator must be bigger than the peak 
voltage and is to be lower than motor back emf [15]. 
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Figure 1. The BLDC motor represented by its equivalent circuit 


The 3-phase voltages for the motor connected to its inverter can be given as follows: 
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where V5, Vps and V,, are the motor stator voltages; R is the stator phase resistance; i,, i, and i, are the stator 
phase currents and e,, e, and e, are the stator back emf’s phase voltage. 
The above phase voltages equations can be extracted in terms of the line voltage equations as [2, 15]: 














Vab R —-R_ 0OTfla L-M M-L 0 4 Ls €a — ep 
Vbc| = | O R —RIfip] + 0 L-M M-L cr Ip] + |&b — &c (2) 
Vea —R 0 R Ic M-L 0 L-—M Ic Eo — €g 




















with considering a balance 3-phase source and 3-phase stator winding, two equations can be extracted due to: 
1,+ip)+i1,=0 (3) 


These two equations are given in terms of phase a and b as follows: 


a: R. Z 1 

aia = — [ola + 3 opp Vab — Cab) + 3G yp oe — be) (4) 
. R. « i 1 

atib = — Foy lb — 37cyp (Va — Cav) + 3G yp oe — &e) (5) 


The back emf's in BLDC motor has a trapezoidal shape depend on rotor position, which is phase- 
shifted by 120 degree. These back emf's equations can be expressed as [2, 15]: 


€a = “: WmF(8¢) (6) 


ep = : WmF(8¢e — =) (7) 
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Ke Tl 
ec = 2 WinF(O¢ + *) (8) 


where k, is the constant of back-emf; w,, is the motor speed and @, is the rotor angle (electrical), as: 
O. = 5 Om (9) 
P: number of the poles and 6,,, is the rotor mechanical angle, as: 
Om = J, Wm dt (10) 
The function F(6,) represents the trapezoid back emf's, for one cycle as: 
10<0®<<= 


1-£(6, -# 26.2 a 


3) 3 
-InsQ<2 (11) 
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F(8.) = 


The torque in the rotor generates power as [17]: 

Da = Gala Fehly F ecle (12) 

Pe = T.Wy (13) 
where T, is the developed torque which can be from equations (6-8, 12 and 13) as follows: 

Te = “t(F(O,)ig + F(Oo — )in + F(Oe + Dic] (14) 

Generally, the motion equation includes: 


dwm 
Te — Ty = JE + kW (15) 


where T;, is the load torque, k, is T, constant, ks is constant for viscous friction and J is the inertia for the rotor. 


3. ROTOR SPEED AND LIMITING LINE CURRENT CONTROL IN BLDC MOTOR DRIVE 

As a comparing among the motor types, the motor of the BLDC motor has more problematic control 
strategy. The well-known methods of controlling a BLDC motor are three: control by hysteresis band, also 
controlled by a variable DC-Link voltage and finally control by PWM for the switching signals [15]. Typically, 
the common control system is a closed dual-loop. The interior loop is responsible for the torque adjustment as 
well the current to the proper value, whilst the other circuit is the outer speed circuit, that’s hired to control the 
speed of the motor. Figure 2 shows the graphical diagram of the drive for the BLDC motor. Subtract the 
reference speed by rotor speed to produce the speed error signal which 1s fed to a speed controller that its result 
is confined by the upper reference current. 
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Figure 2. Graphical diagram of the drive for BLDC motor PWM control system of the speed 
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A reference current value is derived from the following reference torque relationship [2]: 


Tre 
leet -_ a (16) 


The current of DC-link equivalent signal is synthesized of the 3- phase stator's currents thereafter 
subtracted of the output of a speed controller. Thereafter, the signal of error is fed directly to a controller of a 
torque. Here the output controller of torque is fed directly to the form of PWM to compare it through a triangular 
waveform as the comparative output is a low signal or a high signal that turns as the chopping signal into the 
inverter as appeared in Figure 4. The speed of BLDC motor immediately alterations by a duty cycle modulation 
of the firing signals inverter switches leans on the error signal of the controller. In this method of control, the 
motor will switch on and off through a high average, the frequency of chopping 1s steady but the error of control 
will change the length of a duty cycle [2]. The actuality of the constant frequency lets the acoustic filtering and 
electromagnetic noise easier. The frequency of switching is usually taking (20-50) kHz where the maximum 
frequency of switching will give a smoother torque as a result of a low distinction in the current. The control 
system suggestion is consists of suggested sensors of currents to apply a requisite 3-phase currents to snaffle a 
current of DC-Link. This is advantageous for staying the current of a motor at the coveted value through the 
starting and during the speed step variation by existence the controller of a torque. By this way, the signals of 
firing are controlled the length of the duty cycle for a PWM; speed controlling is obtaining. 


4. HARMONY SEARCH OPTIMIZATION METHOD 

Harmony search (HS) is a phenomenon-mimicking metaheuristic algorithm which has been proposed 
in [16]. HS is a novel approach inspired by the musical operation of searching an ideal case of harmony as in 
music improvisation where the pitch of each musical instrument is improved to accomplish this better state of 
harmony. The pitch determines the quality of music, just as the fitness function value determines the quality of 
solutions. Furthermore, the improvisations of musician are analogous to native and global seek sketch in 
optimization strategies. HS has been used to solve various types of optimization problems [17]. It is showing 
significant improvements over other heuristics. A lot of modified HS algorithms have been studied to enhance 
the performance of the original version [18]. 


4.1. Technique of Optimization 

The Techniques of optimization are broadly classified into deterministic techniques of optimization 
and heuristic techniques of optimization. In this paper, it has implemented the Harmony Search Algorithm 
(HSA) as a technique of optimization because it’s an efficient and effective inspired algorithm from nature. 
Executive steps of that HSA have been shown in Figure 3. HSA was suggested by Zong Wo Geem et al., [16] 
in 2001. Seeking for better and new harmony it's called improvisation of harmony. Process of improvisation 
accomplishes in a band of musical where any musician has several tools to play. Any band from the musician 
can be seen as a solution space decision variable. Every tool has several pitch zones which depicts decision 
variables range [16]. The musical band particular iteration, at a set time through the process of improvisation, 
leads to the musical harmony that is afterward the Harmony Memory (HM) is updated. HM is the structure of 
a core data utilized in HSA and it has best solution matrix vectors. HM size describes vectors processed counts 
with each other. The parameters of HSA are selected manually which are shown in Table 1 of this paper for an 
optimal solution of the PID parameters. 
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Figure 3. Procedural steps of HSA 
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4.2. HSA Executive Steps 

The HS algorithm has five executive steps viz: process initialization, initiation of harmony memory, 
improvisation of unprecedented harmony, HM updating and test a stop criterion. All these proceedings are 
briefly qualified below [19]: 

a) Process of initialization: 

Define the lower and upper limit of PID variables x? < x; < x¥, where xi is any gain of the PID 
parameters, L and U is the lower and upper limit respectively of these parameters. Harmony Memory Size 
(HMS) value decision, Harmony Memory Consideration Rate (HMCR) value decision, Pitch Adjustment Rate 
(PAR) value decision, Fret Width (Bandwidth) (FW), Fret Width Damp Ratio (FWdamp) and Define maximum 
iterations 

b) Initiation of the HM: HM matrix is generated randomly for HMS out of the following equation. 


where x(i,j) Shows initial HM jth element. Rand () function produces random values amidst 0 and | whereas 
lj and Uj are lower and upper bound respectively. 

c) New harmony process of improvisation: a process of Improvisation for the new harmony is started by 
step 3 of HSA. This process consists of three sub-steps: HMCR, random selection and PAR. generate a 
number randomly between 0 and 1 with the help of the rand (). Insomuch as the randomly produced value 
1S minimal than HM consideration rate thereafter the value is chosen of that particular position. Insomuch 
as arandom number is maximal than HMCR thereafter generate a new random number. 

This process will appear in the equation below: 


x(randj) if randb( )is < HMCR 
= l, + rand().U; —1, else ct) 
j sae as 
The elements of memory are selected from the HMCR step, are further amended with PAR. The PAR 
equation has appeared in the below. 


(19) 


‘ rand( ). bwj if rand( )is < PAR 
a 


v) else 


d) HM update: when new process of harmony improvisation completes then updates a new harmony in the 
HM worst among already present HMs is removed from HM. 

e) 5) Checking the criteria of stopping: when improvisation process arrives at its iterations maximum 
number or halting criteria of Integrate Time Absolute Error (ITAE) of the multi-objective function as 
shown in equation (20) below, then calculation process of HSA cut out. Otherwise, the process repeats 
from step 3 and step 4. 


F= ft le dt+ ft» le|dt (20) 


where e, is the speed error and e; is the current error. 


5. PROPOSED SIMULINK MODEL FOR BLDC MOTOR BASED ON HAS TECHNIQUE 

The set of equations (1-15) are used to be modeled by the software of Matlab/Simulink (R2015a) at 
Bogaki-Shampine (ode3) with the solver of the fixed step kind by primary sample time is 3.125*10° seconds 
to perform BLDC motor supplied from a 3-phase inverter. This model 1s shown in Figure 4. The BLDC motor 
used in this work has the specifications given in Table 2. 

The proposed Simulink model for BLDC motor including the dual loops is shown in Figure 5. 
According to the schematic diagram given in Figure 2, each control loop includes a PID controller. The 
parameters of these controllers are optimized by the HS method. The outer loop controls the motor speed while 
the inner loop controls and limits the line current to a value that remains the motor in safety operation condition. 
The objective function given in equation 20 is represented in Simulink as shown in Figure 6. The optimized 
parameters for PID controllers are given in Table 3. 
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Table 2. BLDC motor specifications Table 3. Optimized parameters of PID 

Type of parameter value Unit Parameters Values 
Output power 1200 W Motor speed controller 
Line current 16 A Proportional (k,) 13.940 
Terminal voltage 76 Vv Integral (k;) 0.1268 
Full load torque 2.9 N.m Derivative (ka) 0.02022 
Current at no load 0.66 A Current controller 
The speed at no load 3500 Rpm Proportional (k,) 9.829 
Motor torque constant 0.207 Nm/A Integral (k;) 6.599 
The resistance between two 0.110 QO 

terminals 

Motor voltage constant 0.207 V.s/rad 

L-M inductance 0.6 mH 

Rotor inertia 1.7x1107? kg.m? 

Motor viscous damping factor 1.3x1074 N.m s/rad 

Poles number 8 
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Figure 4. BLDC motor controlled by 3-phase inverter Simulink model 
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Figure 5. Proposed dual control loops for BLDC motor Simulink model 
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Figure 6. Multi-objective function block diagram 


4. SIMULATION AND RESULTS 
The simulation results show that the PID speed and current controllers’ gains are shown in Table 3 


which obtained from HSA during run the system with its parameters shown in Table 1. 

Figures (7-9) shows the speed response, stator current for phase (a) and developed torque respectively 
for the BLDC motor drive system with and without current controller by starting with 2000 rpm reference 
speed and applying 2.9 N.m at 0.5 second during | second simulation time. The performance of the BLDC 
motor by changing the reference speed by trapezoidal input and starting directly with a full load is shown in 
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Figure 10 fore speed response and Figures (11 and 12) for developed torque and current respectively surely 
with the current controller. These results ensure the validity of the model, the speed, current and torque reach 
to its new values quickly. 
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Figure 7. Speed response with and without a current controller 
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Figure 8. The current phase (a) with and without a current controller 
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Figure 9. Developed torque with and without a current controller 
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Figure 10. Speed response by trapezoidal reference speed. 
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Figure 11. Developed torque response by trapezoidal reference speed 
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Figure 12. Current response by trapezoidal reference speed 


5. COMPARISON BETWEEN CONTROLLING METHODS 

The comparison is done concerning to the operation and results from the output waveform of the 
performing models as shown in Table 4. The proposed controlling method has a complex hardware structure 
but it gives a better load rejection. And the 2" current loop will dissipate small power through the CT's of the 
currents measurement. 


Table 4. Comparison between controlling methods 


Feature Without current controller With current controller 
Rise time Faster Slower 

Starting and running response Good Better 

Ripple High Low 

Accuracy Less More 

Investment cost Small Higher because of the 24 loop 


6. CONCLUSION 

The current controller was setting the current and overcome the high starting current. In this paper, 
the PID parameters were obtained using HSA it's found is easier and faster to get the result than other used 
optimized technique. This is used to control the BLDC motor speed and current also can provide better 
performance of the motor variables. The benefit of the using HSA 1s to reach a better response of the speed, 
overcoming the load rejection and good robustness system during trapezoidal reference speed change at direct 
starting load. The results of different BLDC motor variables are plotted, discussed and compared. The current 
controller has more rising time in speed but good load rejection and less oscillation. The result of current and 
torque shows the more ripple in the ramp change of the speed but this disadvantage can overcome by re-running 
the optimization algorithm in this state. 
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